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Abstract 

We study the generalized matrix model which corresponds to the n-point toric Virasoro conformal 
block. This describes four-dimensional M = 2 SU(2) n gauge theory with circular quiver diagram by 
the AGT relation. We first verify that the generalized matrix model is obtained from the perturbative 
calculation of the Liouville correlation function. We then derive the Seiberg-Witten curve for AT = 2 
gauge theory as a spectral curve of the generalized matrix model. 
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1 Introduction 

Recently, an interesting conjecture has been proposed in [T] that there are two equivalent ways for 
describing two M5-branes wrapped on a Riemann surface. One of the ways is a four-dimensional TV = 2 
superconformal SU (2) quiver gauge theory, which is constructed depending on the number of genus and 
punctures of the Riemann surface [21 [3j. The other is the Liouville theory defined on that Riemann 
surface. The proposal, which is often called AGT conjecture pQ, is that the Nekrasov partition function 
[1] of such quiver gauge theory can be reproduced by the n-point conformal block of the Virasoro algebra, 
where n is the number of punctures. This conjecture has been generalized to the higher rank gauge group 
|H1 [7J IB], non-conformal case [TU1 E], and to the case in the presence with the surface and loop 
operators [H]-[2I]. ( See also [23 [H] for the M-theory considerations.) A proof has been given in [2"T] 
and [28] for SU{2) gauge theory with adjoint and (Nf — 0, 1,2) fundamental hypermultiplets by using 
the recursion relation [2l?1 1301 13TI 152] . 

It was discussed in [33] that the AGT conjecture for the case on a sphere can be understood through a 
matrix model. The matrix model expression of the conformal block on a sphere is given by reinterpreting 
the Dotsenko-Fateev integral representation of it as the (beta-deformed) matrix integral [Ml [33] with a 
logarithmic potential. When the central charge c — 1+6Q 2 (Q — b+l/b) equals one, or equivalently, b = i, 
the conformal block is represented by the usual matrix model. The matrix model technique, in particular 
large N limit, is useful to show that the gauge theory result can be reproduced [351 133 [Ml I3S] • (See also 
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[30] for 1/N correction.) It has been discussed that direct integral calculation leads to the instanton (q- 
)expansion of the Nekrasov partition function (and the corresponding expansion of the conformal block) 

m M MM M MM- 

While the understanding of the AGT conjecture for a sphere through the matrix model has been 
developed, a counterpart for a Riemann surface with higher genus is not yet fully understood. In this 
paper, we study the following integral: 

Z = f\{ dz, Y[ 0*(zi - Z] )- 2b2 cxp \ J2 2mk lo & e *( z j ~ w k) + 4nriazA j , (1.1) 

J i=i i<j \ gs j=i \k=i J J 

where the function #*(z) is defined in terms of the Jacobi theta function 9\{z\t) and the Dedekind eta 

function 77(1") as 8*(z) = q~ 1 / 12 8i(z\T)/rj(T) with the fixed modulus q — e 27rlT of the torus. This integral 

is an extension of the beta-deformed matrix model corresponding to the conformal block on the sphere 

to that on the torus with n punctures, where z — w is replaced by 9*(z — w) 33 . Q A closely related free 

field representation for toric conformal block is discussed in [15] . 

Although this integral cannot be written in terms of a usual matrix integral, it can be seen as a 

"generalized matrix model" , whose "eigenvalues" Zi live on the torus. Rigorously speaking, the integral 

is defined on the cover of the torus because the integrand itself is not completely doubly periodic. However, 

the positions of cuts, at which the eigenvalues are placed, are indeed doubly periodic. When we regard 

the integral variables z% as the eigenvalues, the product of theta functions Yli<j @*{ z i ~ z j)~ 2b can be 

regarded as a counterpart of the Vandcrmonde determinant. This form is expected from the propagator of 

the two-dimensional free theory on a torus. The parameters are the momenta of the vertex operators 

while Wk arc their insertion points. The parameter a and the filling fractions are identified as the internal 

momenta in the conformal block. The parameters in the integral are related as 



E 



m k + bg s N = Q (1.2) 



due to the momentum conservation law. 

The AGT conjecture indicates that the n-point conformal block on a torus is related to the Af = 2 
superconformal SU(2) n gauge theory with a circular quiver diagram. The parameters rrik are identified 



tin 1421 . a slightly different integral was studied. Instead of introducing a factor iiriazj as in 1 11. 111 . ./Vi integrals over A- 
cycle and N2 integrals over B-cycle are introduced in order to reproduce the expected number of free parameters. However, 
it was also discussed in }421 that their results do not completely reproduce the conformal block although it is very close. 
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with the masses of the bifundamental hypermultiplets and the internal momenta of the conformal block 
correspond to the Coulomb moduli parameters. We see that such relation can be partially understood 
through the integral (jl.ll) in the large N limit: On one hand, we show that the integral is derived from 
the perturbative calculation of the correlation function of the Liouville theory. On the other hand, we 
confirm that the spectral curve of the integral can be identified with the Seiberg-Witten curve of the 
quiver gauge theory. 

The organization of this paper is as follows. In section [5J we review the M-theory construction of 
the four-dimensional superconformal SU(2) n gauge theory and see the M-theory curve which describes 
the low energy effective theory of the gauge theory. In section [31 we see that the integral representation 
(|1.1[) is derived from the Liouville n-point correlation function on a torus by perturbative calculation. 
We then discuss that the proposed integral (11.11) corresponds to the toric conformal block. In section 4, 
we consider the large N spectral curve of the integral and identify it with the curve obtained in section 
O We conclude with discussions in section [5] 

2 M-theory curve of SU(2) quiver gauge theory 

In this paper, we consider TV = 2 superconformal SU(2) n gauge theory with circular quiver. This type 
of quiver gauge theory can be constructed as a worldvolume theory of D4-branes suspended between 
NS5-branes in type IIA string theory [2], as depicted in figure [Tj The D4-branes occupy the x ' 12 ' 3 and 
x directions and the NS5-branes occupy x ' 1 ' 2 ' 3 and x 4 - 5 . The x 4 ' 5 coordinates are combined into the 
complex one v = x 4 + ix 5 . The difference of the ^-coordinates of D4-branes in neighboring intervals of the 
NS5-branes is identified with the mass parameter of the bifundamental hypermultiplet. The x 6 direction 
is compactified and is periodic. 

The curve which describes the low energy effective theory of this quiver gauge theory can be obtained 
by the M-theory up-lift by adding periodic a; 10 direction and considering the hypersurface in v,z = 
x 6 + ix 10 space. Since both the x 6 and a; 10 directions are compactified, these represent a torus E. We 
denote this torus by Weierstrass form 

y 2 = Ax 3 - g 2 x - g 3l (2.1) 
where these coordinates are related to z in terms of Weierstrass elliptic function as x = p(z) and y = p'(z). 
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Figure 1: Brane configuration of the SU(2) 4 gauge theory with circular quiver. The D4-branes are 
suspended between the NS5-branes. 

The curve in X = C x E which the M5-brane wraps on is identified with the Seiberg-Witten curve. For 
the SU(2) n case, the form of it is F(x, y, v) = with 



The positions of n NS5-branes are translated in M-theory to the points in the torus E in which fi and /2 
have simple poles [2]. The residues of f± are interpreted as the mass parameters rrik of the bifundamentals. 
(When we discuss poles or residues, we use the local coordinate z = x 6 + ix la .) 

Note that the above consideration corresponds to the case where the total sum of the hypermultiplct 
masses vanishes: 



In order to include the most generic case with non-zero sum m of the masses, we have to consider non- 
trivial bundle X m — > E such that x e — > x e + 2irR, v — > v + m, rather than the trivial bundle X. By 
this choice, the constraints on j\ and f% are as follows: let Pk € E (k = 1, . . . ,n) be the positions of 
the singularities coming from the NS5-branes. Then, /i and fa have order 1 and 2 poles respectively at, 
say pi, and have simple poles at P2, ■ ■ - ,Pri- (Note that the sum of the residues of f± should vanish and 
therefore the residue at pi of /i is — Y^k=2 m k-) Furthermore, the double pole of fa at z\ disappears 
when we use a good coordinate v = v + S ^ . 

However, as pointed out in [3], the problem can be simplified by eliminating the linear term in v and 
writing the curve as follows 



F{x,y,v) = v 2 - fi(x,y)v + fa(x,y). 



(2.2) 



n 




(2.3) 



k=l 




(2.4) 
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where we have shifted as v — > v + # and therefore 02 = 4 



ji. It follows from the above constraints 



that 4>2 has double poles at Pi,P2> • • • >Pn an d the coefficients are simply the mass parameters m 2 ,. The 
Coulomb branch parameters are included in less singular terms: we can introduce n — 1 parameters Ck by 
adding the functions which have simple poles at pk because the sum of the Ck should vanish. We are also 
free to add the constant term in 4>2- These correspond to n Coulomb moduli parameters. Indeed, the 
most singular part of <p2 can be represented by p(z) and the less singular part is by elliptic £ function: 
p = — j-C- Therefore, the curve is v 2 — <p2 with 



where Zk are the values of the local coordinate z = x 6 + ix w at the points p^- This curve can be seen as 
a double cover of the torus with n punctures whose positions are z = z^. The Seiberg-Witten one form 
is represented by vdz. 

3 Prom Liouville theory to generalized matrix model 

In p] , it has been proposed that the n-point conformal block on a torus can be identified with the Nekrasov 
partition function of Af = 2 quiver gauge theory discussed in the previous section. This relation for a 
torus was further studied in [51] [27]. The identification of the parameters of the two is as follows: the 
momenta of the vertex operators, the internal momenta and the complex structure moduli correspond, 
respectively, to the mass parameters of the hypermultiplets, the Coulomb moduli and the gauge coupling 
constants in the gauge theory. The Liouville correlation function can be obtained by integrating the 
contribution from the three-point function and the conformal block over the internal momenta. In this 
section, we explain how the integral representation (jl.ip appears from the full correlation function, based 
on the discussion of [52[ [33] . 

The n-point function of the Liouville theory on a torus is formally given by the following path integral 



n 



n 




(2.5) 



\k=l 




) 



Liouville on T 2 



V<t>{z,z)e- S[4,] Y[e 



k=l 



,2m k (p(wk,Wk) 



(3.1) 



where the Liouville action is given by 




(3.2) 
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under the flat background metric. We choose the insertion points Wk such that they satisfy 

n 

^rrikWk = 0, (3.3) 
fe=i 

which does not break generality due to the translational invariance of the torus. 

We divide the Liouville field into the zero mode and the fluctuation (f>(z,z) — 4>o + (j>(z,z). By 
integrating over ^>o, we obtain 

. v->n / n \ „ / „ \ — y\ • • 

Z-, k = l b I . rvi . \ ' / / •, — 

_ 1 "6 



.4 



26 



u-=i 



l r ( £ ^ ] / ^(z,z) e - So ^ ( / d 2 z e 2b ^ 



fe=i 



2m k cf>(w k ,w k ) 



,(3.4) 



where So is the free scalar field action. When 

n 



(3.5) 



the correlator diverges due to the factor T(—N). Up to this divergent factor, the n-point correlation 
function A is given by a perturbation from the free theory 



, JV r N I N 



n 

fc=i 



a 2m k (f>(w k ,w k ) 



(3.6) 



free on T 2 



The condition Q3.5P ensures the momentum conservation in the free theory. 
The £-point function of the free theory on a torus is given bjj] 



n 



Jki4>(zi,Zi) 



iC* 2 (t)(2tt)S(Y^ h) J] "M^i { Zij \t) exp 



free on T 2 



T2 



(3.7) 



where = Zi—Zj, r is the modulus of the torus, T2 is its imaginary part, and C^ 2 = (4-7r 2 T2) 2 |?7(t)| 2 ■ It 
factorizes into holomorphic and anti-holomorphic parts by introducing an additional integral as [54} l55ll53] 



n 



i=1 I free on T 2 

/>oo 

2i\T 1 {T)\- 2 8{Y^h) / da 



n 



q a exp I ~2iri'^^kj Z ja 



, (3.8) 



where q = exp(27rir). We recover (|3.7p by explicitly carrying out the Gaussian integral over a. This 
expression is more suitable for associating the Liouville n-point function to the holomorphic generalized 
matrix model. 



§See the review |53| and the references therein. 



Here, we introduce 

0.(2) = 



77 (r) aj - 



9 m )(l 



(3.9) 



where g = exp(27ri-r) is fixed, for later convenience. Using the explicit expression (|3.8|) for (|3.6|) . we find 
that the n-point function A of the Liouvillc theory reduces to the following integral 



A = C{r,m k ,b) Yl 



(Wkl 



exp 



l<fc</<n 

N n 



da \q\ 



-2a 



N 



T- 



i=l 



N 



-26^ m k log 9* (zi - wk) - 2b 2 log 6* (z%j ) - Amba z% 

l<i<j<N 



i=l fc=l 



i=l 



,(3.10) 



where ro. 



Wj , and we have introduced the factor in front of the z integral as 



C(r,m fc ,6) 



S(0)\v(r)\- 2 \q- 1/24 v(rT 4 ^ mk2 ~ 4b2N . 



(3.11) 



The discussion above is valid even for finite N. However, it is not straightforward to factorize the 
integrals over the torus into holomorphic and anti-holomorphic integrals for generic N. In order to 
proceed to the next step, we evaluate the integral (|3.10l) in the large N limit. From the momentum 
conservation (13.51) . we see that m k = 0(N) as b = 0(1). Although the integral variable a runs over the 
whole range of imaginary numbers, dependence on a of the exponent appears when a = O(N). We see 
that all the three terms in the exponent in (I3.10P are all 0(N 2 ). In the large N limit, the integral (|3 . 10[) 
is evaluated at the critical points of the exponent of the integrand. The conditions for the criticality of 
the exponent are given by 



> rn k — r+fr> „ ; \ + 2ma = 

^ H ¥ {Zi-W h ) ^ 



t{ z ij) 



(3.12) 



k=l ~»> jjii -'^31 

where 9'^z) = d z 6*(z). The conditions obtained from the Zj-derivatives are just the complex conjugate 
of (|3.12[) . It is remarkable that the conditions for criticality are separated into holomorphic and anti- 
holomorphic equations, which indicates that the integrals over the torus in (|3.10p can be factorized into 
holomorphic and anti-holomorphic integrals in the large N limit. 

If we do not have the second term in (|3.12p . we have n possible critical points for each variable Zi, 
assuming that the parameters nik are generic. We expect that the n critical points are "diffused" to form 
line segments due to the second term, similarly to the case of the usual large N matrix model. Then, the 



solutions of (13.121) are labelled by the filling fraction v k = bg s N k , in which Nk out of N variables Zi take 
the value on the fc-th line segment. Here, we have introduced the parameter g s , where g s N is finite in 
the large N limit, for later convenience. 

We define the contributions from the solution labelled by {vk} to the holomorphic integral as 

N 



Z(q,w kl m k ,a,v k ) = j Y\_ dz i \\ 



-2b A 



l<i<j<N 

exp ( -— ( 2m k log9*(zi - w k ) + Amazi | ) , (3.13) 



where we have rescaled the parameters as m k — > m k /g s and a — > a/g s . The paths of the integrals are 
defined such that only the solution of (|3.12j) labelled by the fixed filling fractions {vk} contributes to 
the integrals. By regarding the factor Yii<i<j<N ^*( z ij) 2b as generalization of the Vandermonde 
determinant, we see that the holomorphic integral p.!3|) is "the generalized matrix model" with the 
action 

n 

W(z) = ^2m k log9*(z - w k ) +4maz. (3.14) 

k=l 

The integral in (|3.10p is then obtained by integrating p.!3|) and its complex conjugate over the filling 
fractions. Thus, in the large N limit, the n-point function A of the Liouville theory can be written as 

/>zoo i* 

A = / da dv\ ■ ■ ■ dv n -\ 
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(q- 1/U v(r)) II 0,(w kl )-^ m '^\ q ~ a2 ^Z( q ,w k ,m k ,a^ k ) 



Kk<l<n 



,(3.15) 



where we have used that C(r, m k , b) is approximated as \q 1 / 24 ?7(t) I 4 ^" fc ^ 9s . 



The total n parameters, a and the independent filling fractions i>i, ■ ■ ■ , v n -i, can be identified with 
the n internal momenta eti, • • • , a n in the conformal block B(q, w k , m k , a k ) |33j . Under this identification, 
we see from p. 151) that Z(q, w k , m k , a, v k ) corresponds to the conformal block as well as the holomorphic 
contribution from the three-point functions. 

In the next section, we relate the generalized matrix model (I3.13|) with the Af = 2 quiver gauge theory 
in the previous section. Before going into it, let us see the remaining parts in (|3.15[) here. As discussed 
in PQ, the gauge coupling constants = exp(27rirj) (£ = !,•■■ ,n) of the SU(2) n quiver gauge theory is 



related to the modulus q — exp(27rir) of the torus and the insertion points Wi of the n-point function of 
the Liouville theory as 

qi = e 2m ^-^\ 92 = e 2 ™^-^), q n _ 1 = e 2 * i( - w ^- w «\ q L q 2 — q n =q. (3.16) 

Using (|3.9[) and the definition for the Dedekind eta function q~ x ^ 2 ^r\(r) = fj" t (l — q m ), the factor in 
front of Z(q, Wk,mk, a, Vk) in (|3.15|) is rewritten in terms of qi as 

-2j2 k 7n k 2 /g 2 s 



l<k<l<n 
n oo 

= n ■ ■ ■ <n-i) mkmi/9 °) n n (! - ^ ■ ■ • *+*r 2mimi+fc+i/ff = (3.1?) 



k<l i=l k=0 

where the subscripts of m, and qi are considered modulo n. This factor has a close form as the U(l) 
factor, which is explicitly written in pQ. Furthermore, the factor q~ a l 9s in (|3.15|) corresponds to the tree 
level contribution to the prepotential. 

In this section, we have derived the generalized matrix model from the Liouville correlation function 
by explicitly carrying out the perturbative calculation. In [33] , it was discussed that the action (|3.14[) of 
the generalized matrix model is also expected from a geometrical argument in topological string theory 
in the context of the AGT conjecture. 

4 Spectral curve of the generalized matrix model 

In this section, we derive the spectral curve of the generalized matrix model, which is introduced in the 
previous section: 

/ 1 \ f N ( b N \ 
exp I— 5-F] = / n^ cx P ^^(Ai) + I>S(0*( A *- A i)r 26 > 

where W(A) is given by (|3.14|l . As stated previously, the paths of the integrals are determined such that 
they realize the given filling fractions. For generic parameters, the action W(A) has n critical points. In 
order to make the discussion as generic as possible, we will calculate it with a generic action with n critical 
points and, at the final stage, substitute its specific form (|3.14p . We will see that the Seiberg-Witten 
curve (|2.4p with (|2.5p appears as the spectral curve of this generalized large N matrix model. 



9 



In the large N limit, the problem of the integration in (I4.1j) reduces to calculation of the critical point 
of its exponent. The prepotential is given by 



4-F=--^T^(A J )-^6 2 log^(A J -A J ), (4.2) 



where each eigenvalue satisfies the condition of criticality 

0=^W'(X i )+2bJ2 9 ^ff^-. (4.3) 

It is natural to assume that the eigenvalues are distributed around the critical points of PT(A) but 
in the form of line segment, similarly to the usual matrix model. We denote these line segments as Ck 
where k = 1, • • • n. We do not assume that Ck are on a real axis. However, we assume that Ck do not 
include the singular points we, at which the action W(A) diverges. Suppose that Nk eigenvalues are on 
the line segment Ck, where Nk satisfies 2fc=i = N. 

Here, we introduce the density of eigenvalues p(A) which has non-zero value only on the line segment 
C\- Outside of these regions, we define that p(X) = 0. The density of eigenvalues is normalized as 
J c dXp(X) = bg s Nk = Vk- Using the variables introduced above, the prepotential and the condition for 
criticality are written as 



T=\ dXp(X)W(X) + dX' dXp(X)p(X') log#*(A — A'), (4.4) 

= W(X)' + 2 / dX'p(X') 6 ; iX ~ X 'l (4.5) 

respectively. 

In order to solve this, we define the resolvent as 

R(z)= [ dXp(X) ^^- . (4.6) 

Since the behavior of the theta function around the zeros is 9*(e) ~ Ce + 0(e 3 ), 01(e) ~ C + 0(e 2 ) 
for small e, the structure of the singularity is similar to that of the usual matrix model if we focus on 
the fundamental region of the torus. The resolvent has cuts at the line segments Ck- Also, the filling 
fractions are obtained by integrating the resolvent along the cuts as 



— I dzR(z). (4.7) 
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The significant difference from the usual hermitian matrix model is that the resolvent has pseudo 
periodicity due to the theta function. Using the identity 

— log 6*{z + m + nr^r) = -2%in + log 0*(z,r), (4.8) 

which holds for arbitrary integer m and n, the pseudo periodicity of the resolvent is shown to be given as 

R(z + m + tit) = R(z) - 2ninbg s N. (4.9) 

We see that the resolvent is completely periodic for the A-cycle, which is parallel to the real axis, but 
some constant is added for the £?-cycle. This pseudo periodicity gives great restriction on the possible 
form of the resolvent below. 

On the line segments Ck, the resolvent is expected to behave as 

R(z + iee l ^ z) ) + R(z-iee lv(z) ) = 2 f dX' p(X') 6 *} = -W'(z), (4.10) 

Jj2 k c k V*(z-A) 

R(z + iee lip{z) ) - R(z - iee lv(z) ) = I d\' p(\') 9 *) Z ~ ^ - -2mp(z), (4.11) 

where we take real number e infinitely small and f(z) is properly defined such that z + iee 1 ^ or 
z — iee 1 ^^ does not go across the cuts Ck when z moves along Ck- The integral in (|4.10|) is principal 
integration, which is given as an average of integral along the path above the singularity and that below 
the singularity. The resolvent should be determined such that (|4.10|1 and (14.11| arc satisfied for z € Ck, 
together with the periodicity (|4.9p . 

A candidate of the solution for (|4.10p is 

Ro(z) - -\w\z). (4.12) 

However, it does not reproduce the correct structure of singularity expressed in (|4.11[) . We need singular 
contributions: 

R(z) = R (z) + R(z) siDg , (4.13) 

where (|4.10p and (|4.11[) impose 

R{z + tee lv ^) sing + R(z - iee lv ^) sing = 0. (4.14) 
R(z + iee l ^ z) ) sing - R(z-iee lv ^) sing = -2mp(z), (4.15) 
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The above discussion is valid for generic action W(z). In the following, we use the concrete form 
(|3.14p to determine the resolvent R(z). Since Rq(z) is given by 



where we used the identity (|4.8j) and the momentum conservation ()3.5j) . (Note that the parameters nik 
and a have been rescaled.) Remarkably, the pseudo periodicity (|4.17l) for Rq(z) exactly agrees with that 
of the resolvent R(z) in (|4.9I) . Therefore, we find that R(z) s i ng is exactly doubly periodic. 

On one hand, we find that the resolvent R(z) has no singularity except for the cuts in the regions Ck 
from the definition (|4.6p . On the other hand, we see from the explicit form (|4.16p that Rq(z) has simple 
poles at z = Wk + m + nr and their residues are —m^. Therefore, -R(z)sing must have both the cuts in 
the regions Ck and simple poles with residues rrik at z — Wk + m + nr. 

From (|4.14|) . we see that the sign of i?(z) s i n g changes when crossing the cut. Thus, R(z) 2 ing has no 
such discontinuity and all the cuts disappear. Since R(z) 2 - lng has double poles at z = Wk + m + nr, 
the Laurent expansion of R(z) 2 lng at z = Wk has not only (z — Wk)~ 2 terms but also (z — Wk) terms 
due to contact terms. We denote the coefficients of such terms as Cfe. Such doubly periodic function 
is uniquely determined, up to adding a constant, as a linear combination of the Weierstrass £ function, 
which has one simple pole in the fundamental region, and of the Weierstrass p function, which has one 
double pole in the fundamental region. This uniqueness is followed from the Liouvillc theorem, which 
ensures that a function holomorphic on the whole complex plane is only the constant function. In other 
words, we cannot add non-trivial function to the doubly periodic function without changing its structure 
of singularity. 

From the discussion above, we have derived that i?(z) s i ng 2 must be of the form 




(4.16) 



fe=i 



we see that its pseudo periodicity is given by 



Ro(z + m + nr) = Rq(z) — 2irinbg s N, 



(4.17) 



R[z) + -W'(z)j = R(z) sing 2 = mip{z -w k ) + J2 c k((z ~ v>k) - C. 



(4.18) 




where the coefficients satisfy the condition 



n 




(4.19) 



k=l 
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so that (|4.18p is doubly periodic. This is the spectral curve of the generalized matrix model. This spectral 
curve (|4.18p coincides with the form of the Seiberg-Witten curve (I2.4p with (|2.5p . 

4.1 One-point function 

In the following, we concentrate on the case of one-point function and determine the constant in the 
spectral curve under the approximation that \a\ is large enough. The spectral curve (|4.18p for n = 1 
reduces to 

R(z)sin S 2 = m 2 p{z) - C, (4.20) 
where we used the condition (|4. 19[) . From (|4. 1 5[) , the density of eigenvalue is written as 

P (X) = (±)!V™ 2 P(A)-C, (4.21) 



where the sign should be determined so that —ip(X)dX is positive real on the integral path C 



i- 

The equation of motion (|4. 5[) for n = 1 reduces to 

Now, suppose that ia 3> 1. In order that the eigenvalues satisfy the equation of motion, at least either 
the first or the third term of (|4.22[) must be large enough to compensate the second term. The first term 
becomes large if an eigenvalue is placed close to the singularity A — 0. The third term becomes large if 
the distribution p(X') takes large value around the pole A' = A of the integrand. Since p(X') is constrained 
by the normalization condition 

/ d\'p(X') = bg s N, (4.23) 
Jd 

p(X') must take large value only around the region close to A' = A and rapidly decrease as it goes apart 
in this case. However, since p{X') is given by (|4.21l) . this occurs only if A is close to the singularity A = 0. 
Thus, in any case, every eigenvalue A must be distributed very close to the singularity A = 0. In other 
words, the cut C\ of the resolvent is placed very close to the singularity z = and its length is very short. 
When A ~ 0, the eigenvalue density (|4.2ip can be approximated as 



^Va^- c2 < (4 ' 24) 
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c 


^ - 


-1/c -e 


£ 1/C 



Figure 2: Path Ci 



Figure 3: Path C[ 



where we put C = m 2 c 2 for future convenience. Therefore, in this approximation, the endpoints ai, 61 
of the cut C\ (p(ai) = p{b\) = 0) are 

1 



ai,bi ~ ± : 



Also, the equation of motion (|4.22l) can be approximated as 



777 

— + 2lTia - 

A 



Ci 



Since the solution of the classical equation of motion is given by 

A 



7777 



(4.25) 



(4.26) 



(4.27) 



it is expected that the cut C\ is placed around this point. Taking account the symmetry, most natural 
possibility is that the endpoints a\, b\ = ±l/c are on the real axis and the path C\ goes beyond the origin 
as in figure [H 

For later calculation, we analytically continue the eigenvalue density p(X), which has non-zero value 
only on the cut, to the whole complex plane and we denote it as p(z). However, if we define p(z) by the 
function in (14.241) under the usual convention that the function f(z) = yfz has cut on the negative real 
axis, it changes the sign across the imaginary axis. Since the original eigenvalue function p(X) does not 
have such discontinuity on Ci, we rewrite (|4.24p and define p(z) as 



p(z) EE — y/1 {CZ)\ 
■KIZ 



(4.28) 



so that it is smooth on the path C\. The sign was determined so that —ip{z) is positive at z ~ im/a. 

At this stage, let us check the normalization condition (|4.23p . By using the analytically continued 
function (|4.28|) , we can change the path of the integral from C\ to C[ as in figure [3] because there are no 
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cuts or singularity between the path C± and C' x . Then, we can explicitly calculate the integral as 

ps p rl/c 

dXp(X) = / dzp(z) + / dzp(z) + / dzp(z) 

Ci J-l/c J z=ee ie ,6:tt^0 Js 

1/c 1 

dz[p{z) + p(— z)] — — x 27riRes z= o p(z) 



(4.29) 



Thus, taking account the momentum conservation bg s N + m = 0, we find that the normalization (|4.23l) 
is satisfied for arbitrary c and imposes no constraints. 

We now use the equation of motion (|4.26l) at A = 1/c, which is the endpoint of the cut G\. 

cm + 2ma+ [ dX' P } X \, = 0. (4.30) 
Jd 1/c - A 

By changing the path of the integral similarly to the previous calculation and by using (|4.28|) . we find 
that the integral in the equations of motion can be explicitly calculated to give 

/ dX \f X '\, =-(l + *)mc. (4.31) 
Thus, we have determined the coefficient c as 

2ira 

c= , 4.32 

m 

which indicates C — (2ira) 2 . Thus, the spectral curve (|4.20p is 

R(z) 2 sing = m 2 p(z) ~ (2™) 2 . (4.33) 

In the limit \a\ S> |m|, it is straightforward to check that 

/ R(z)sin g = 2nim } / R(z) sing = 2nia, (4.34) 

where j3 is the small cycle around the origin while 7 is the ^4-cycle of the torus. These are the expected 
property, which the Seiberg-Witten curve should satisfy. Thus, we have shown that the spectral curve 
of the generalized matrix model coincides with the Seiberg-Witten curve including the constant term in 
this limit. 
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5 Conclusion and discussion 



In this paper, we have studied the generalized matrix model, which is a Dotsenko-Fateev type integral 
representation of the toric conformal block. This generalized matrix model is shown to be naturally 
derived from the perturbative calculation of the n-point function of the Liouville theory on the torus. 

We have shown that the Seiberg-Witten curve of the corresponding gauge theory is derived from 
the generalized matrix model as a spectral curve. We have confirmed that the constants of the spectral 
curve of the generalized matrix model with n = 1 agrees with those of the Seiberg-Witten curve for 
TV = 2* theory in a large internal momentum limit. Our results suggest that the AGT relation between 
the Liouville theory on a torus and the four dimensional M = 2 quiver gauge theory can be understood 
through the generalized matrix model. 

One of the future directions is to determine the undetermined constants in the spectral curve in more 
generically. Then, we go on to the next step to calculate the prepotential and compare the results with 
the conformal block or with the Nekrasov partition function. 

Other direction is to consider a finite N correction. We have seen the large N limit has played an 
essential role to derive the generalized matrix model from the Liouville n-point function. Thus, such an 
extension seems to be quite non-trivial. In a sphere case, it is known that the matrix expression for the 
conformal block is valid even for finite N [40] - [47]. So, it would be quite interesting to study whether 
such an expression is possible also for the torus case. In particular, the derivation of the loop equation 
would strongly help to calculate finite N partition function [SHI E3 EE] ■ 

The generalized matrix model which we have considered in this paper can be seen as an elliptic 
extension of the Selberg integral. The properties of this have been studied in [SHI |HS]. It would be 
interesting to study the generalized matrix model from this direction. 
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